Superfluidity versus Disorder in the Discrete Nonlinear Schrodinger Equation 
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We study the discrete nonlinear Schrodinger equation (DNLS) in an annular geometry with on- 
site defects. The dynamics of a traveling plane- wave maps onto an effective "non-rigid pendulum" 
Hamiltonian. The different regimes include the complete reflection and refocusing of the initial 
wave, solitonic structures, and a superfluid state. In the superfluid regime, which occurs above 
a critical value of nonlinearity, a plane-wave travels coherently through the randomly distributed 
defects. This superfluidity criterion for the DNLS is analogous to (yet very different from) the 
Landau superfluidity criteria in translationally invariant systems. Experimental implications for the 
physics of Bose-Einstein condensate gases trapped in optical potentials and of arrays of optical fibers 
are discussed. 
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Studying the interplay between nonlinearity and disor- 
der has become a fundamental issue of the last decades in 
the study of many physical and biological systems (both 
discrete and continuous) |Q. It is well known that non- 
linearity or disorder may lead to localized excitations - 
'solitonic' structures || and Anderson localization' j|, 
respectively. However, the dynamical properties of the 
system, when both nonlinearity and disorder are present, 
are still challenging theoretical investigations. This prob- 
lem is, moreover, of central experimental relevance, since 
impurities can be reduced but never completely elimi- 
nated. In particular, it can be asked if random defects 
will (and, if so, how) destroy the propagation of traveling 
plane waves or localized excitations (allowed by the non- 
linearity), and what are the conditions for crossing from 
a "superfluid" regime with propagation (and coherence) 
preserved (due to, for instance, a large nonlinearity), to 
a "normal" regime with disorder induced damping. 

Here, we consider the dynamical properties of the 
DNLS in an annular geometry and in the presence of im- 
purities. We choose the DNLS for two main reasons: 1) 
it has all the required ingredients: nonlinearity, disorder 
and discreteness; 2) real physical systems, such as opti- 
cal fibers and Bose-Einstein condensates in deep optical 
lattices, map onto the DNLS and provide an ideal exper- 
imental framework. The annular geometry is paradig- 
matic for studying the nature of superfluidity [Q, and 
allows a clear comparison between the discrete and con- 
tinuous (translationally invariant) limits. 

With regard to optical fibers, a typical experimental 
setup is an array of one dimensional nonlinear coupled 
waveguides [||. As the light propagates through the ar- 
ray, the coupling induces an exchange of power among 
the single waveguides. In the low power limit (i.e. when 
the nonlinearity is negligible), the optical field spreads 
over the whole array. Upon increasing the power, the 
output field narrows until it is localized in a few waveg- 
uides, and discrete solitons can finally be observed H||. 



The evolution of E n (r), the electrical field in the nth 
waveguide, as a function of the position, r, is governed 
by the DNLS Eq. ([j]) . In this case A is proportional to the 
Kerr nonlinearity and the on-site potentials e n are the ef- 
fective refraction indices of the individual waveguides. In 
|Q a linearly growing e n was realized, which allowed the 
observation of Bloch oscillations. 

Another significant physical system described by the 
DNLS is a Bose-Einstein condensate (BEC) gas con- 
fined in a deep optical lattice. A one dimensional optical 
lattice can be created experimentally by a far-detuned, 
retro- reflected laser beam p|-|TT|. The condensate wave- 
function vb(r, t) obeys the Gross-Pitaevskii equation |L^ | 
which is a (continuous) nonlinear Schrodinger equation, 
with the nonlinearity arising from the interatomic inter- 
action. When the heights of the intra-well barriers of the 
periodic optical potential are much higher than the con- 
densate chemical potential, the system can be mapped 
onto the DNLS Eq.(Q) (T3), with ip n the condensate am- 
plitude in the nth well and e„ proportional to any exter- 
nal field superimposed on the lattice. In || a coherent 
output of matter waves was created by a vertical optical 
array (with the gravity gradient providing e„ cx n). In 
pOf , the magnetic trap and the laser beams were turned 
on in a superimposed harmonic magnetic field (e„ oc n 2 ), 
allowing the direct observation of coherent (Josephson) 
condensate oscillations governed by Eq. ([j]). 

The DNLS is (in dimensionless units): 



.dip n 



■~(lpn-l + tpn+l) + (e n + A\lp n | 2 )i/>„, (1) 



where A is the nonlinear coefficient and n = 1, • • • , N 
(N number of sites). In the physical systems we have 
discussed, the defects e n can be spatially localized or 
extended. For instance, the impurities in optical fibers 
can be induced by different (possibly random) effective 
refraction indexes of the guides or with varying spatial 
separations between them. In BEC's the defects can be 
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created with additional lasers and/or magnetic fields; the 
presence of a thermal component can also be phenomeno- 
logically modeled, in some limits, by a random distribu- 
tion of defects. 

We consider, first, the DNLS with a single impurity 
e n = e 5 n ,n at the site n, and study the propagation of a 
plane wave ip n (r = 0) = e lkn . In the following we assume 
A > (which corresponds to a repulsive interatomic in- 
teraction in BEC's, as is the case for 87 Kb atoms). Note, 
however, that Eq.(|]) is invariant with respect to the 
transformation A — > —A, e n — > — e„, and ip n — > jfj n e lvn . 
Since we consider periodic boundary conditions (due to 
the annular geometry), we have k = 2nl/N with I integer 
(Z = 0,1,---,7V-1). 

In the translationally invariant limit of the DNLS, 
given by the continuum nonlinear Schrodingcr equation 
(CNLS), a well know argument suggested by Landau im- 
plies that superfluidity occurs when the speed is smaller 
than the sound velocity (for weak perturbations). A sim- 
ple derivation of the Landau critical velocity in the CNLS 
was recently proposed in B considering an annular ge- 
ometry with a single (small) impurity. The key point was 
to map the problem of the propagation of a plane wave to 
a Josephson-like Hamiltonian. The superfluid regime is 
allowed by the nonlinearity, which provided an effective 
energy barrier against the creation of elementary exci- 
tations with momentum k + q,k — q (with q arbitrarily 
small) which would dissipate the energy of the incident 
wave having momentum k. 

This scenario is completely changed by discreteness. 
First, it is well known that when cosfc < the system 
becomes modulationally unstable jL4j. Stability analysis 
reveals that the eigenfrequencies of the linear modes be- 
come imaginary driving an exponential growth of small 
perturbations. This modulation instability disappear, for 
A > 0, in the CNLS limit. Let us consider, then, the case 
in which cos k > 0. In the absence of the impurity, super- 
positions of rotational states with opposite wave-vectors 
k, — k are degenerate. The defect splits the degeneracy by 
coupling the two k, —k waves, very much as the tunneling 
barrier does in a double well potential, with "left" and 
"right" localized states. Therefore, the relative popula- 
tion of the two waves oscillates according to an effective 
(generalized |H|| ) Josephson Hamiltonian. These Joseph- 
son regimes are preserved as far as the splitting in energy 
induced by the defect ~ e is much smaller than the en- 
ergy gap between different rotational states ~ 27rsinfc. 
In this limit, we can write the wavefunction ip n {' r ) as 

i) n (r) = A(T)e lkn + B( T )e- lkn . (2) 

In the following, we set A, B = ^n A ^ B (T)e l<t ' A - B( - T \ 
z = riA —tlb and <f> = <pA — 4>b- We will compare the 
numerical solution of (Q) with the analytical solution of 
(||) obtained from the ansatz (||). 

The two-mode Eq. (||) can be extended to the case of a 
time-dependent, arbitrary (including random) distribu- 
tion of defects, with e replaced, as shown below, by an 



effective impurity strength. Furthermore, when the ini- 
tial wave function is given by the sum of multiple waves 
4>n{0) = j A? £tkj ™ ' ^ ne ansatz (||) can be straightfor- 
wardly generalized so long as the quasi-momentum dis- 
tributions peaked around kj do not overlap. The collision 
of a soliton with a single impurity has been studied, from 
a different perspective, in [ fl5| |. A numerical analysis of 
the propagation of plane waves across a segment with 
defects was made in 

Let us now derive the equations of motion. We define 
an effective Lagrangian as L = J2 n tynipn ~ ^ where 

n = EJ-i(« + i + rjM+i) + e n i I 2 +4 1 ^ I 4 ] 

(both the Hamiltonian Ti and the norm | ip n \ 2 ~ N 
are conserved). The Euler-Lagrange equations ^J^- = 
for the variational parameters qi(r) — ita.b, (Pa,b 
give the following equations: 

z = \/l — z 2 sin <h (3a) 

N y ' 

^ = iw^ C0S ^ + Az < ( 3b ) 

with the replacement <f>+2kn — » <\>. The total (conserved) 
energy is: 

^=^-|v / T^cos^ (4) 

and the equations of motion (||) can be written in the 
Hamiltonian form z = — 0^ and 4> = with z and (j) 
canonically conjugate variables. 

The Eqs.(||) have been studied in very different con- 
texts, including polaron dynamics, where the dimcr 
Eqs.(||) has been solved analytically jTjj, and in the 
Josephson dynamics of two weakly coupled Bosc-Einstein 
condensates |lq] . Eqs.(||) are those of a non-rigid pendu- 
lum: <p is the angular position and z its conjugate mo- 
mentum. The non-rigidity of the pendulum is due to its 
momentum dependent length. 

The pendulum phase portrait, z~4>, has been studied 
in flisfl . Let us briefly recall the main results. We have a) 
oscillations around < (j> >— and < z >— (0-states); 
b) oscillations around < z >^ with running phase 

< (f> >oc t (self-trapped states); c) oscillations around 

< z >= and < cf> >— n (n-states); d) oscillations about 

< z >^ and < 4> >= it (self-trapped ir-states). Here 

< • • • > stands for a time average. The stationary points 
of Eqs.(|), z = 0, 4> = (i.e. A = B) and z = 0, <$> = it 
(A = —B), correspond to time-independent solutions of 
Eq.(|l]), ip n — 2 cos kn and tp n = 2i sin kn, respectively. 

To understand the meaning of these regimes in our sys- 
tem, we observe that the angular momentum is L(t) — 
'EB(^C+r c - c -) = 2Nzsink : < z >= implies that 
the wave is completely reflected, and < z(t) » (or 

< z(t) >< 0) that the wave is only partially reflected by 
the impurity. The latter regime is given by a complete 
rotation of the pendulum about its center, and can be 
considered as a self-trapping of the angular momentum. 
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Equivalently, there is an effective energy barrier which 
forbids the complete reflection of the incident wave, and 
preserve its coherence. The observation of a persistent 
current is associated to a superfluid regime of the DNLS 
equation. 

The situation changes in the (quasi-)continuum limit 
(which should not be confused with the N — > oo limit). 
In this case phonons can be emitted only with quasi- 
momentum close to k, a condition which allows the 
applications of the Landau superfluidity criteria. The 
crossover between continuum and discrete limits is a very 
interesting one: We emphasize that the two-mode dy- 
namics is crucially related to discreteness and nonlincar- 
ity, and disappears in the CNLS limit. Yet, there is a 
striking analogy: in both cases (in the Landau and in the 
" pendulum" criterion) the phonon emission out of the in- 
cident wave (which, therefore, dissipates its energy) can 
be inhibited by an effective energy barrier. The key dif- 
ference lies on the corresponding spectrum of the emitted 
phonons, which leads to a completely different dynamics. 

From the effective Hamiltonian Q we find a critical 
value A c for the pendulum oscillations about its center 
given by 

H(4>(0),z(0)) = 2e/N: (5) 

when A < A c , z oscillates around 0. When A = A c , 
asymptotically z(r) — > and with A > A c , < z(r) >^ 0. 
In Fig.l we plot the average value of the normalized angu- 
lar momentum L(t)/Lq for different values of A/A c and 
z(Q) = 1, 0(0) = 0. The numerical solutions of Eq.ffl) 
arc in agreement with the two- mode approximation (0), 
dashed line. In the inset of Fig.l we plot the normalized 
angular momentum vs. time for different A/A c . 

The fixed points of Eqs.(^) can be found by solving 
for i = 0, = 0. In particular, we have the non- 
trivial stationary (solitonic) solutions = (2m + 1)tt, 

z = ± yl — (with m integer). In Figs.2(a)-2(d) we 

compare the numerical and pendulum solutions for the 
normalized angular momentum and the phase, in the 
cases of tt oscillations (a-b) and it stationary points (c-d). 

The previous discussion can be extended to the case of 
many impurities: replacing the ansatz Eq.(Q) in Eq.(Q) 

we obtain: H = 4^- - j|Vl — z 2 J2 n e« cos (0 + 2kn). 
After some algebra, the effective Hamiltonian becomes: 

. Az^ 26 / ;r , , , 

H = — --Vl-2 2 cos(0 + a) (6) 

with e and a given by the Fourier transform of the defects 
distribution: 

ee ,a =^e n e Mn (7) 

n 

The critical value A c is given by Eq.(||) with the re- 
placement e — > e; 0o — > 0o + a (19) ■ It is also clear, 
from Eq.(R), that the system becomes transparent for 



some particular distribution of defects. For instance, 
with an extended, step-like barrier (e„ = constant for 
n\ < n < T12) of length 10 sites and with 2k = it/5, we 
have e = 0. 

All the predicted regimes discussed so far have been 
found to be in agreement with our full numerical analysis. 
We considered different uniform random distributions of 
defects (e.g. all the e„ positive or negative or with zero 
mean value) . The critical values of the nonlinearity found 
from the numerical solution of Eq.(Q) and the compari- 
son with the theoretical prediction for A c from Eq.(||) is 
shown in Fig. 3. In the inset of Fig. 3 we plot L(t) /Lq as a 
function of time for various A and a random distribution 
of defects e n . The excellent agreement between numer- 
ics and the solution of Eq.(||), and the robustness of the 
two-mode ansatz in the presence of an arbitrary distri- 
bution of defects, opens to the possibility of studying 
the competition between disorder (and Anderson local- 
ization) and nonlinearity from a new perspective. Eq.(|^) 
is analytically solvable, yet still it contains all the essen- 
tial ingredients to investigate the details of the superfluid 
- normal transition in the DNLS with impurities. 

To conclude we briefly discuss the limits for recover- 
ing the CNLS equation (in an annular geometry) from 
the DNLS Eq.(g). Writing A = 2mg Q /h 2 N, e n = 
V n mL 2 /TiN 2 and t = mL 2 /TiN 2 T, with V n = V(x = x n ) 
the defect potential in x n , L the length of the annulus and 
r the dimensionless time entering in Eq. ([!]), the CNLS is 
obtained in the limit N — > 00. In particular, the crit- 
ical value for the pendulum oscillations Eq.(^) becomes 
A c = V n mL 2 /TlN 3 — > 0. Therefore, approaching the con- 
tinuous limit, the DNLS pendulum regime collapses to a 
(strongly) self-trapped state. This prevents the emission 
of phonons with opposite momenta with respect to the 
incident wave, whose energy will be eventually dissipated 
on a much longer time scale, according to the Landau ar- 
gument. 
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FIG. 2. Normalized angular momentum (a) and phase (b) 
vs. time for a 7r-state (A/A c = 0.5). In (c) and (d) we plot the 
same quantities for A/ = (2e/N) /a/1 — z 2 (0) (the stationary 
solution of Eq.(|) with </> = tt) and for A/A/ = 0.95, 1.10. 
In all cases solid (dashed) lines are for the numerical (varia- 
tional) solution of Eq.([j]) with e = 0.01, N = 100, z(0) = 0.5, 

(f>(0) = TV. 
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FIG. 1. Average value of the angular momentum L(t) 
(normalized to the initial value Lo) vs. the nonlinear co- 
efficient A/Ac, (A c = ie/N), with e = 0.01, N = 100, 
z(0) = 1. The filled circles are the numerical solutions of 
Eq. (|lj) , the dashed line is obtained from equations (||). Inset: 
normalized angular momentum vs. time for different values 
of A/A c = 0.5,0.75,1,1.5,25, respectively, corresponding to 
(o),-.,(e). 0(0) = 0. 




FIG. 3. Critical value A c vs. e for different random dis- 
tributions of defects (and initial values z(0) — 1, <f>(0) — 0). 
Black circles: numerical solution of Eq. (jlj) . Solid line: the an- 
alytical prediction A c = 4e/N with e given by Eq.(Q). Inset: 
angular momentum vs. time with a random distribution of 
defects for different values of A/A c = 0.45, 0.90, 1.01, 10, 1000 
(corresponding to (a), • ■ ■ , (e)). The sum of the strengths of 
the random impurities is 0.1. 
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